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ABSTRACT 

We study how massive halos respond to perturbations. Through numerical 
solution of the coupled linearized Boltzmann-Poisson equations, we find that halos can 
transmit and amplify disturbances over large distances within galaxies. In particular, 
as Weinberg has noted, the halo provides an excellent medium for transmitting 
disturbances from the outer regions of galaxies to the inner regions. The dipolar 
response typically dominates and, in some cases, is very weakly damped. The monopole 
and quadrupole response can also make significant contributions. Overall, the results 
support the notion that disk structure can be excited by transmission of noise from the 
distant halo. The strong dipolar response suggests that halos play a role in producing 
lopsided disks. These results also suggest that, during formation from initial collapse, 
halos do not settle rapidly into static equilbrium; instead, they may continue to ring 
and interact with other embryonic components of proto-galaxies. 

Subject headings: stellar dynamics - galaxies: individual (Milky Way) - galaxies: 
haloes - galaxies: kinematics and dynamics 

1. Introduction 

Galaxies typically inhabit noisy neighborhoods. Noise- or irregularites in the gravitational 
field- usually appears in the form of large satellites in the outskirts of galaxies (e.g. Zaritsky 
& White 1994). It is also expected in the form of continuing infall and inhomogeneity in the 
dark matter distribution. Yet most of this noise appears to lie at rather large distances. We 
might suppose, therefore, that it does not disturb the tranquility of the inner regions of galaxies, 
particularly galactic disks. 

However, massive halos dominate most galaxies. They extend to large distances and, 
in some sense, provide a medium which potentially fosters communication between spatially 
separated regions, as Lynden-Bell (1985) first pointed out. Recent work by Weinberg (1995,1998) 
elucidates the dynamical behavior of this connective medium: the halo can respond to forcing by 
a satellite such as the LMC and efficiently transmit disturbances into the inner regions of galaxies. 
Consequently the apparent spatial isolation of inner regions is balanced by the halo's ability to 
transmit disturbances: the centers of galaxies seem to hear the noise, loud and clear. 
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Given the variety and ubiquity of disk structure, this might not come as a surprise. 
Galactic disks prominently display warps, lopsidedness, bars, spiral arms and other asymmetries. 
Theoretical work on disk structure shows that individual features are short-lived and must 
be continuously re-excited. Neverthess, it has been difficult to pinpoint the culprit since 
many candidate disturbers such as the LMC lie too far from the disk to excite structure 
through direct tidal forcing. Halo transmission and amplification seems to provide at least one 
necessary mechanism for continuously exciting disk structure (see Nelson & Tremaine 1995 for a 
comprehensive review of disk dynamics). 

Unfortunately, it is difficult to study the response of halos to moderate amplitude 
perturbations. Ideally we could simply perform simulations. However, current N-body techniques 
suffer from resolution limitations due to particle number, drift in orbit integration and, in the case 
of expansion codes, inability to conserve momentum. These problems suppress linear collective 
effects which we expect to dominate the response in this regime. 

An alternative approach, known as the 'matrix method', was originally developed to calculate 
dispersion relations of disk systems by Kalnajs (1977) and has been adapted by Weinberg (1989) 
to study spherical systems. The original implementation is best suited for stability analysis and 
determining the growth rates of unstable modes (e.g. Weinberg 1991; Saha 1992) because the 
dispersion relation is defined in the upper-half of the complex frequency plane. Recently, Weinberg 
(1994) modified the method to identify modes in stable systems using elegant numerical techniques 
to continue the dispersion relation into the lower half-plane. 

In the present work, we further modify the matrix method in order to calculate the evolution 
of a perturbed system in time. In essence, we solve the coupled linearized Boltzmann-Poisson 
equations as an initial value problem to determine the linear response of the model system. Our 
approach is well suited to the study of stable galaxy response because the time-dependence and 
relative amplitudes of excited modes are uniquely determined for arbitrary external perturbations. 
Weinberg (1989,1995,1998) has taken a similar approach in calculating the response in the 
asymptotic limit of a periodic perturbation. 

Herein we present the results of calculations using this method. Our results support the 
findings of Weinberg (1994,1995,1998) and demonstrate that halos can amplify perturbations 
and absorb energy into weakly damped oscillatory modes. Perhaps as importantly, the results 
demonstrate the need for high precision numerical calculations to study the subtle effects which 
determine the structure and evolution of galaxies. The outline is as folllows. We first derive the 
modification of the matrix method in §[| We then discuss the numerical implementation of the 
method and the battery of tests used to verify its accuracy in §||. The first section on results, 
§[|, details the response of King models to point-like perturbations. Although not presented in a 
direct astrophysical context, these results demonstrate the variety of effects arising in the models 
and provide a convenient framework for categorizing the physical behavior. In the second section 
on results, §[|, we present an example of a fly- by to illustrate how effects arise in an astrophysical 
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context. Several possibilities for further investigation are discussed in §^j. 



2. Formalism 
2.1. Hamiltonian 

We choose coordinates centered always on the unperturbed (monopole) potential. This 
introduces an indirect term into the Hamiltonian due to the acceleration of the reference frame. 
The resultant Hamiltonian then only describes the internal reaction of the unperturbed system to 
the perturbation. For tidal encounters, the indirect potential directly cancels the I = 1 component 
of the perturber's gravitational potential. However, for interpenetrating encounters, the term 
appears explicitly. 

The full Hamiltonian is 

H = i|v| 2 + $ (|r|) + $ s (r,t) + <S> e (r,t) + ^. r ^H + AH. (1) 

The indirect term is given by the linear potential dV c /dt ■ r, where dV c /dt is the mean acceleration 
of the unperturbed system. Defining the total potential perturbation $>i = <£ s + 3> e as the 
respective sum of the response and direct external potentials, we can write 

dV c If If, d$ ft) . 

-ir = -wj drp ^ = wj drpi ^ (2) 

where p\ = p s + p e is the respective sum of response and external source densities. The indirect 
term introduces a small modification of the matrix equation for the £ = 1 component of the 
response. In these coordinates, the response contains no barycentric shift. 



2.2. Linear theory and matrix method 

The matrix method for determining the stability of stellar systems was originally developed by 
Kalnajs (1977). Our derivation is related to Weinberg's (1989) adaptation of the matrix method 
to spherical systems. 

We consider the effect of weak perturbations on an equilibrium stellar system with distribution 
function F(l) (hereafter DF), where the vector I denotes the actions associated with unperturbed 
orbits in the potential. Jeans' theorem states that the equilibrium DF does not depend on the 
angle variables w (Binney & Tremaine 1987). The evolution of the perturbed DF /(I,w) is 
described by the linearized Boltzmann equation 

df dH df OF dAH 

~dt + ~~&T ' ~ ~cW ' ~d^T ~ ' ( ' 
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where AH is the perturbing Hamiltonian defined above. The response potential $ s is related to / 
through Poisson's equation 

V 2 $ s = 4vrG J fdv = 4irGp s . (4) 

In the non-self-gravitating approximation, one sets $1 = <E> e . 

The perturbed quantities, / and AH, may be expanded in Fourier series in the action-angle 
variables of the unperturbed system (Goldstein 1980; Tremaine & Weinberg 1984). Following 
their convention, we define I\ as the radial action, I2 as the total angular momentum and 
I3 as the z-component of the angular momentum. Motion in the unperturbed Hamiltonian 
Hq(1) = ^|v[ 2 + <I>o(l r l) is given by I = constant, w = wt + wo. In particular, the Fourier expansion 
of the perturbation becomes 

AH{v,t) = ^A# k (I,t)exp(ik- w), (5) 

k 

where 

A# k (I,t) = ^3 j dwAF(r,t)exp(-ik-w). (6) 

Tremaine & Weinberg (1984) give expressions for the action-angle series of a multipole potential 
expansion: we introduce these formulae below. 

Substituting into equation @ yields the relation 

— — — +ik-uf k = ik-u—AH k (I,t), (7) 

since dH^/dl = u> defines the vector of stellar orbital frequencies and dF/dl = udF/dE for an 
isotropic DF. We can construct the inhomogeneous solution to equation (fjj) from the homogeneous 
solution (right-hand side zero); therefore 

dF C l 

/ k (I, i) = tk • u— / dt' exp[ik • u (f - t)]AH u (I, t'). (8) 
at/ J -00 

Summing over k and integrating over velocities yields a formal expression for the perturbed 
density: 

r dF f l 

p s = / dv^exp(fk • w)(ik • w)— / cft'exp[ik • u(t' - t)]AH k (l,t'). (9) 

J , dF J -co 

k 

To determine the response density and potential, we first introduce biorthonormal expansions 
(e.g. Kalnajs 1977; Weinberg 1989) for the density, 

p s = Y,aT (t)tf m (r)Y em (Sl), (10) 

£mj 

and potential, 

& = J^aj" 1 (t)u' m (r)Y em (n)- (11) 

Imj 
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The perturbing potential $ e has an analogous expansion in terms of coefficients denoted bj m (t). 
The condition of biorthonormality relates the basis functions through Poisson's equation 

V 2 uf™ = 4nG4 m , (12) 

and imposes orthogonality, 

"drr 2 uf m di m = 5ij. (13) 



4ttG 

To substitute into equation @, we first rewrite these expansions in action-angle variables 
using expressions from Tremaine & Weinberg (1984). The coefficient 



A# k = 5>f 3 (i)+&f(i)] 



Wg (I) + ^- P « m X k (I)S ei 



VgktkM, ( 14 ) 



where cos/? = J3/I2 and V^ 2 /% 3 (/3) is defined in terms of rotation matrices. The second term in 
brackets on the right-hand side arises from the indirect potential. The Fourier coefficients have 
the definitions 

1 Z* 71 " 

W^(I) = — dwi ex.p(-ikiwi)u e j k3 exp[i/c 2 (V ; — ^2)], (15) 



and 



1 f n 

X k (I) = — / d^iexpf— i&iiiMrexp — 102)]. (16) 
The angle ip — W2 has the definition 

dr(^ 2 - I 2 /r 2 ) 

ifj-w 2 = (17) 
^p(i) ^2[£7-*(r)] -/|/r 2 



Lastly, 



^ = Jdrr*d^. (18) 



Now substituting these expansions into equation (||), taking the inner product with respect 
to ui m *Y e * m (ft), where 



im*v* 



Y? m {n) = VU^Wti 1 ) exp(-ik' • w), (19) 



and noting that drdv = dldw, we obtain a set of coupled integral equations for the expansion 
coefficients: 

«' m W = -^cJ dIdw ^ E (*-o,)|drex P [ik.u,(r-t)][af (r) + 6f (r)] 



^kH(/3)^ 2 fc 3 (/3)^k"(I) 



<(I) + ^fX k (I)% 



exp[i(k - k') • w] (20) 
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Integrating over angles yields (27r) 3 5j c k'- To integrate over (3, we employ the orthogonality relation 
of rotation matrices (Edmonds 1960): 



2 

dcos pVff k2ks {(3)V £k2k3 (P) = - \Y ek2 (±n,0) | 2 See = C ek2 5 e , £ , 



where 



Ceko 



2 2k 2 -l (^_fc 2 )l 

vr 2 (£ + k 2 )\ 





r[i(i + fc 2 + i)] 

T[\{l-k 2 ) + l] 



I + ki even 
, £ + k 2 odd. 



(21) 



(22) 



Changing variables from I\ to E, where dl\ = dE/co\, we define the matrix kernel 

4tt 



ij 



(2vr) 3 f JdJdF ^ ... ... . a 

' dE — ^ C/fe, (zk ■ w) expfik • w(r - t)]W£ 



loi dE 



. (23) 



This reduces equation (pC|) to the form 

af m (*) = E T *-i^(T - t)[af» (r) + 6f (r)]. 
The equations can be solved iteratively (Tricomi 1957). 
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To derive the matrix equation in the frequency domain (Weinberg 1989), we take the two-sided 
Laplace transform of equation (24), so that 



a(s) =M(a)[a(s) + b(a)], 



where the matrix 



(2-)! 
4vrG 



dE 



JdJ dF 
oj\ dE 



ik • 



, S + IK ■ LV 



4ir 



(25) 



(26) 



2.3. Green's function and linear response operator 

Define the Green's function for the density response in the following way: 

p s (r,t) = J dr'dt'G p (v,v'-t-t') P e (r',t'). (27) 

G p depends only on t — t' because the system is translationally invariant in time. The notation 
G p indicates that the Green's function determines the response density since we can define an 
analogous relation for the response potential. Introducing a point source density 

p e (r',t')=S(r'-r p )5(t'-t p ), (28) 
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yields an equation for the Green's function 

PiUM) =G p (r,v p ;t-t p ), (29) 

where the notation p rp ,t p indicates that the response is due to a point source at a particular 
location. 

The linear response operator has a definition similar to that of G p (Nelson & Tremaine 1997): 

p s (r, t) = J dr'dt'R(r, r'; t - t')<5> e {r', t'). (30) 

Comparing equations (^?j) and ( |30|) shows that the Green's function and linear response operator 
have the following relation: 

R{r, r; t-t') = G p (r, r; t - t')(4vrG)^ 1 V 2 r , (31) 

The method defined above provides a numerical definition of the Green's function or linear 
response operator for a system. 



3. Numerical implementation and tests 

We use standard numerical techniques to implement the method. There are three basic 
components: the action-angle expansions, the basis function expansions and the solution of the 
coupled integral equations. Angles and frequencies are computed from the orbit equations on a 
120 x 20 grid in E and J, respectively, with k\ < 4. Calculations with k\ < 8 produced the same 
results: higher-order resonance influence damping at small spatial scales rather than the large 
scales of interest here. We use the Hernquist basis for the biorthogonal expansion and truncate the 
series after 30 radial terms. The coupled integral equations are solved iteratively to a tolerance of 
10~ 6 using Simpson's rule on a 1024-point grid to evaluate the lag time quadrature in equation 



(24) for r > t - 200. 



In practice, we first compute the angles, frequencies and Fourier coefficients of the basis 
functions for a particular model and t and dump the results to a file using an unformatted output 
stream. Taking this data as input, we then compute the response kernel in equation ( p4] ) and 
again save the results using an unformatted output stream. Computing the kernel takes the most 
time- roughly 3 hours on a fully loaded SGI origin. However, given the kernel, we can calculate 
the response to any perturbation in short order. 

We performed three tests of the implementation. The simplest was to reproduce the analytic 
adiabatic fluid response calculation for both the isothermal sphere (Murali & Tremaine 1997) 
and a Wo = 5 King model. This provides a test of the spatial response which is independent 
of the temporal response. It ensures that the basis function expansion works properly and 
that a sufficient number of terms have been retained to correctly reproduce the potential and 
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density response. We then tested the accuracy of the spatial and temporal solution by Laplace 
transforming the response coefficients along a strip in the complex plane and plugging these values 
into the dispersion relation, equation (^)Q. The transforms of the dominant coefficients typically 
satisfy the dispersion relation to within at least 10%. 

Finally, as a straightforward but important check, we compare in Figure 1 the results 
of our linear calculation with an N-body simulation performed with the self-consistent-field 
code discussed most recently by Hernquist & Ostriker (1992). The code employs the identical 
multipole-biorthonormal decomposition of the potential used in the linear calculation developed 
above. Thus we may directly compare the time evolution of expansion coefficients computed using 
the two methods. We find good agreement for £ = 2. 

Given the results below, we are particularly interested in comparing the two calculations for 
t = 1. However, we were unable to reproduce the calculated response using an SCF simulation 
with 10 6 particles. As Weinberg (1994) notes, expansion codes do not conserve momentum, so they 
will be hard-pressed to reproduce an off-center response since the center wanders. In particular, 
to defeat wandering, the standard centering algorithm calls for resetting the coordinate center 
to the center-of-mass of the system at each potential recomputation. However, to recognize and 
reproduce a dipole response we should ideally remain in a reference frame which is attached to the 
center of the unperturbed potential (which is no longer at the center-of-mass because of the dipole 
response). Nevertheless, experiments without the centering algorithm also failed to reproduce the 
effect. 



4. Response to point-like perturbations 

In the spirit of calculating the Green's function, we subject the model systems to point-like 
perturbations: 

$ e (r)=5(t)d> p ([r-r p |) (32) 

where & p (r p ) denotes the potential of a point-mass of mass m p located at r p (for definiteness, the 
angular coordinates 9 P = <f> p = 0). Point-like perturbations also approximate Poisson fluctuations 
in equilibrium systems. 

We impose perturbations of harmonic order i = 0, 1, 2 and corresponding m values on a range 
of King models at radii enclosing different mass fractions of the unperturbed system. Table 1 
outlines the computational grid. The notation r x denotes the radius enclosing the fraction x of 
the total mass of the system. The quantity rt denotes the tidal radius or radius enclosing 100% of 
the mass. 

Units are defined such that the total mass M = 1, G = 1 and the total energy E = —1/8. 



Suggested by John Dubinski 
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Fig. 1. — The self-gravitating response of a Wq = 5 King model to an I = 2, m = perturbation 
which is turned on and off exponentially. The four lowest order expansion coefficients are plotted 
in time for the linear calculation (dashed) and N-body simulation (solid) with 10 6 particles. 
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This implies that, at r 1 / 2 , the period of a circular orbit, P 1 / 2 ~ 18. Reaonable physical values 
can be obtained by rescaling the models to halos of mass M = 1O 12 M and limiting radius 

Tmax = 200 kpc. 

We set the perturbation amplitude by adjusting m p . Our choices are arbitrary and serve more 
to illustrate the physics of the response. For distant perturbations M(r p ) > 0.75, we typically 
choose m p = 0.1 while for inner perturbations M(r p ) < 0.75, we typically choose m p = 0.01. Since 
the calculation is linear, we may rescale the results for other choices of m p . 

4.1. Transmission and amplification of disturbances 

We find that significant transmission and amplification of disturbances can occur in the model 
calculations. Since the response to each £ is independent, we consider each separately. 

4-1.1. Monopole response 

The monopole component of the perturbation acts like a spherical shell of mass m p which 
exerts a force only on material at r > r p , the shell radius. This generates a density enhancement 
at r > r p which propagates inward to smaller radii. Figure 2 shows this effect in a Wo = 3 King 
model. We find that the amplitude of the propagating response can be significant well within the 
radius of the perturbation. Since the monopole component of the perturber's direct gravitational 
field exerts no force on material with r < r p , the amplification is infinite. 

Deeper within the Wq = 3 system, smaller perturbations give rise to similar effects in the 
core. For perturbations at the 25% mass fraction, we find a roughly a 1% central response for 
m p = 0.01. The amplitude of the peak central density response decreases with concentration due 
to the increasing depth of the central potential well. 



Table 1: Characteristics of King models and computational grid 
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2.6 
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3.2 


14 
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Fig. 2. — Transmission of spherical disturbance in Wq = 3 King model initially excited at r v = 2.4 
(75% mass radius) for m p = 0.1. Panels show density, potential and respective contrasts as a 
function of radius as indicated. In each panel, solid line shows inwardly growing response at 
t = 1.1 following perturbation. Dotted line shows peak central response at t = 5.8. Dashed line 
shows subsequent minimum at t = 12.8 as disturbance propagates outwards. The perturbation 
subsequently damps away. The peak central density response is roughly 1%. 
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4-1-2. Dipole response 

The dipole component of the perturbation acts somewhat differently from the monopole in 
that it exerts a force at all radii. Some component of this force contributes to the mean motion of 
the system; recall that the indirect potential removes this component in the present calculation. 
Nevertheless, the residual (differential) acceleration excites a strong dipolar response in all the 
models studied here. 

For each of the King models, Table |2] shows the basic response characteristics for perturbations 
at the 75% mass radius. The response travels deep into the system, producing peak density 
contrast at roughly 0.1r p (radius denoted r pea & enclosing M(r pea k) of the total mass). In each 
case, the fractional density response is several %. The response improves with concentration. Also 
note that there may be a substantial delay between the application of the perturbation and the 
peak response in the system's core at t pea k- 

Figure 3 provides a close look at a cross-section of the response in the Wq = 7 system. The 
density at small radii responds strongly to the perturbation. To better understand this effect, we 
show the cross section of the mass perturbation r 2 p s along r p in Figure 4. The perturbation draws 
mass from all surrounding directions. However, mass drawn from inner regions produces a large 
change in the density because of the smaller volume. 

To gauge the importance of self-gravity, it is useful to compute the amplification factor or 
Love number of the perturbation: \ = + $ e ')/<3? e/ , where <3? e/ includes the indirect termP). 
Figure 5 shows \ for the perturbation depicted in the previous figure. The Wo = 7 system achieves 
dipole amplification factors of 10 or so which may persist for long durations in weakly damped 
modes. Such strong amplification will lead to non-linear effects that are not described by our 
present treatment. The amplification factors decrease with concentration: for Wq = 5, \\\ ~ 4; for 
W = 3, |x| ~ 2. 

Similar effects result from disturbances at smaller radius. These lead to stronger effects in 
the core. Table ||| shows the peak central response arising from a point perturbation at the radius 
enclosing 50% of the mass. Here disturbances are transmitted to small radius. Disturbances at 



2 There is apparently some ambiguity in defining y for i = 1: we have chosen "1> S as the total potential perturbation 
and $ e/ as the external potential which produces $ s . 

Table 2: Peak central density contrast for M{r p ) = 0.75 and m p = 0.1 
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Tpeak 




M{r peak ) 


tpeak 
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0.03 


8.1 
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Fig. 3. — Transmission of dipole disturbance in Wo = 7 King model initially excited at r p = 3.2 
(75% mass radius) for m p = 0.1. Panels are as in previous figure with radial profile along r p . 
Solid line shows growing response at t = 3.4 following perturbation. Dotted line shows peak 
central response at t = 8.1. Dashed line shows large subsequent peak after next half-cycle of mode 
oscillation at t = 38.6, indicating the presence of weakly damped modes with very slow pattern 
speeds. The peak inner density response is roughly 5%, dropping to ~ 2% after the next half-cycle. 
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Fig. 4. — The change in mass along the r p axis corresponding to the perturbation depicted in the 
previous figure. At t = 3.4, mass is attracted towards the perturber's original position (r p = 3.2; 
solid). This draws material from surrounding regions, but affects the density small radius because 
of the smaller volume. At t = 8.1, the response peaks (dotted). Finally, after the next half-cycle 
(t = 38.6), the response peaks again (dashed). 
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Fig. 5. — Potential amplification factor x along r p for disturbance shown in figure 

MM (log 



0.5). Solid line shows amplification at peak response (t = 8.1). Dashed line shows amplification at 
peak after next half-cycle (t = 38.6). On the near side of the disturbance, de-amplification of the 
potential perturbation weakens the interior halo potential while, on the far side, de-amplification 
strengthens. The amplification magnitude is roughly 10 and can persist for long times after the 
disturbance disappears, due to the excitation of weakly damped modes. 
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even smaller radius can excite the dominant response directly. 

4-1.3. Quadrupole response 

In these experiments, quadrupole disturbances do not appear to propagate from large radius 
to small radius in the same manner as monopole and quadrupole disturbances. Figure 6 shows 
an example of the response to an £ = 2, m = 2 perturbation in a Wq = 5 King model. The 
main feature is the local density response which then decays away relatively rapidly. This tends 
to produce moderate amplifications of ~ 30% at radii roughly 0.3r p . Overall, the quadrupole 
response appears to be most important in the case of pure resonant forcing (Weinberg 1998). 

5. A realistic example of dipole response 

The above analysis details the physical behavior of the response. To provide a more realistic 
example focusing in the t = 1 response, we take a point-mass perturber with m p = 0.01 on a 
hyperbolic trajectory through a Wo = 7 halo. The perturber initially has velocity v = 0.5 at 
infinity: this is roughly the circular speed at the half-mass radius. Closest approach is at r = 2.6, 
the radius enclosing 68% of the total mass. The maximum velocity is roughly v = 1.0 so that the 
effect of dynamical friction will be small. Figure 7 shows the strong 1=1 response induced by the 
relatively small perturber. 

The encounter excites a very weakly damped mode. Figure 8 shows the time-dependence of 
the induced potential energy perturbation which provides a good signature of mode oscillation. 
For comparison, Weinberg (1994) found an e- folding 2ir/ Im(uj) = 250 and an oscillation period 
2ir/ Re(uj) = 98 in these units for Wq = 7. The response depicted here shows an oscillation period 
P ~ 100 in good agreement. The damping rate is harder to judge because the apparent amplitude 
consists of a complex superposition of interfering modes; nevertheless, the damping rate is low and 
appears to have roughly the same magnitude. Also recall that the precise character of the dipole 
modes will differ because of the difference in Hamiltonian. 



Table 3: Peak central density contrast for M(r p ) = 0.50 and m p = 0.01 
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Fig. 6. — Transmission of quadrupole response in Wq = 5 King model with r p = 1.0 (25% mass 
radius) for m p = 0.01. Panels are as in above figures with radial profile along r p . Solid line shows 
peak response at t = 1.1, immediately after perturbation. Dotted line and dashed lines show 
subsequent decay at t = 3.4 and t = 5.7, respectively. 
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Fig. 7. — The dipole response induced by a point-mass perturber with m f 



0.01 on a 



hyperbolic trajectory in the X-Y plane indicated by dotted curve in each panel. The solid 
square shows the perturber's position when in the field of view. Solid (dotted) contours show 
overdensity (underdensity) in unit logarithmic intervals, log// = —1,-2,-3, about the peak 
at each time. Time interval At = 2.44 separates each panel. The peak density response is 
p s max = 0.022,0.032,0.039,0.042 at r p = 0.2 from left to right, top to bottom. The absolute 
maximum occurs in the final panel after the perturber has flown by and is roughly a 10% 
perturbation to the background density. The mode continues to librate through the center long 
after the perturber has passed. 
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Fig. 8. — Potential energy perturbation as a function of time for Wq = 7 system in Figure ||. The 
potential energy perturbation is only ~ 3 x 10 -4 the potential energy of the entire system. The 
total initial disturbance decays rapidly but a fair amount of the input power persists in a weakly 
damped mode. 
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6. Discussion and Conclusions 

Overall, our results indicate that halos can respond quite actively to perturbations. The 
response appears strong over a wide range of profiles and can efficiently transmit disturbances from 
the outer halo into the inner halo. In particular, the dipole response dominates and is strongest 
for the highest concentration system. The weak damping at high concentration may reflect 
kinematic decoupling of core and halo. The dipole response is still significant at low concentration 
but does not appear to persist in weakly damped modes. At low concentration, the monopole 
contribution can be fairly strong and long-range. The quadrupole response is roughly independent 
of concentration and seems to be more strongly excited through resonant forcing (Weinberg 1998). 

Our results complement recent work by Weinberg (1994). The results appear to be consistent: 
mode shapes are similar; both calculations give roughly the same central displacement due to the 
dipole response; and the weakly damped dipole modes in the Wq = 7 system appear to have nearly 
the same period and similar damping rates in both calculations. Furthermore, the present results 
suggest that it may be difficult to excite the weakly damped dipole modes in Wq = 3 and Wq = 5 
King models discussed by Weinberg (1994). 

More generally, the results suggest that noisy halos play a role in exciting disk structure and 
can drive evolution in the inner parts of galaxies. As an example, recent surveys of disk galaxies 
show a high frequency of lopsided galactic disks (~ 30%; Zaritsky <fe Rix 1997). The strong 
1 = 1 response found in our calculations suggests that the halo has influence in these systems. 
In particular, a weakly damped i = 1 mode with a long period can induce an adiabatic m = 1 
distortion in a disk. Scaling the Wq = 7 results to a 200 kpc halo, we find that the location of 
peak response is roughly 4 kpc, which is consistent with typically observed distortion locations. 
Persistent modes are attractive also because, in many cases, there is no obvious companion to 
the galaxy in question. Thus a previous encounter can excite a long-lived mode that continues to 
slowly distort the disk long after the perturber has gone. Ultimately the disk response depends on 
the relative mass in the mode. This aspect will require further investigation. 

Another significant possibility is that active halos strongly influence the formation of galactic 
disks. The standard scenario is that disks settle adiabatically into static dark halos (e.g. Dubinski 
1994). However, the proto-galactic environment is likely to be extremely noisy, particularly in 
outer regions, so that the halo may perturb the disk by transmitting numerous disturbances into 
the inner galaxy. Moreover, the inner halo may continue to oscillate as it settles after initial 
collapse. Halo oscillations can easily perturb the disk through the time-dependent gravitational 
potential. Conversely, the structural integrity of observed disks set limits on the degree of 
disequilibrium in the proto-galactic halo. 

The above calculations have proved extremely valuable because of their high precision: they 
provide further evidence that the very subtle influence of the halo plays a major role in determining 
the observed properties of disk galaxies. Nevertheless, we emphasize the need to improve and 
perhaps develop alternatives to the methods used herein. The main drawback of the above 
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method is the reliance on numerous expansions and the extreme care required in implementation. 
Moreover, we are limited to systems with spherical geometry whereas galaxies most likely come in 
a variety of shapes. However, commonly used N-body methods still seem far from achieving the 
sensitivity required to reproduce such subtle effects. To reproduce off-center modes might require 
direct N-body codes with at least 10 6 particles- still 1 to 2 orders of magnitude beyond current 
capabilities. 

This research was supported by NSERC and the Fund for Astrophysical Research. I thank 
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Phil Arras, Scott Tremaine, Ira Wasserman, Martin Weinberg, Dick Bond and Peter Goldreich for 
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